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Abstract. We collected photometrical and dynamical data for 73 low red-shift (z<0.2) Radio Galaxies (LzRG) 
in order to study their Fundamental Plane (FP). For 22 sources we also present new velocity dispersion data, 
that complement the photometric data given in our previous study of LzRG (Govoni et. al. 2000a). It is found 
that the FP of LzRG is similar to the one defined by non-active elliptical galaxies, with LzRG representing the 
brightest end of the population of early type galaxies. Since the FP mainly reflects the virial equilibrium condition, 
our result implies that the global properties of early-type galaxies (deflning the FP) are not influenced by the 
presence of gas accretion in the central black hole. This is fully in agreement with the recent results in black hole 
demography, showing that virtually all luminous spheroidal galaxies host a massive black hole and therefore may 
potentially become active. We conflrm and extend to giant ellipticals the systematic increase of the mass-to-light 
ratio with galaxy luminosity. 
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1. Introduction 

It is well known that the global properties of early- 
type galaxies are fairly well described through a three 
dimensional space of observables which, besides the ef- 
fective radius re and the corresponding average surface 
brightness <fie>, involves the central velocity dispersion 
(Tc (Djorgovski & Davis 1987, Dressier et. al. 1987). In 
fact, ellipticals and lenticulars are found to lie in a sur- 
prisingly tight, linear region of this space (the so called 
Fundamental Plane; hereafter FP) that has been shown to 
be close to the plane defining the virial equilibrium con- 
dition over the assumption of a rigorous homology among 
galaxies (Faber et. al. 1989). 

The systematic deviation of the FP from the one de- 
fined by the virial equilibrium, has been interpreted as 
due to changes of the Mass to Light ratio (M/L) (due to 
different stellar content or to differences in dark matter 
distribution; Pahre et. al. 1995, Mobasher et. al. 1999), 
as well as to the breakdown of the homology assumption 
(e.g. kinematic anisotropy; Caon et. al. 1993, Graham et. 
al. 1996, Ciotti et. al. 1996, Busarello et. al. 1997). In this 
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context it is therefore important to determine whether the 
FP has universal validity among ellipticals. In particular it 
is not yet clear if active and non-active galaxies do follow 
the same FP, as it would be expected in view of the re- 
cent results on Black Holes (BH) demography in galaxies 
(Gebhardt et. al. 2000, Ferrarese & Merritt 2000). 

A relevant case is that of low redshift extragalactic ra- 
dio sources that are associated with massive early-type 
galaxies. A number of optical studies have been carried 
out in the last decade to investigate the relationship be- 
tween radio emission properties and the characteristics of 
the host galaxies (see e.g. Govoni et. al. 2000a and ref- 
erences therein). Comparison of the optical properties of 
radio galaxies with those of normal (non radio) ellipticals 
have shown that the probability to exhibit radio emis- 
sion increases with the optical luminosity of the galaxies 
(Scarpa & Urry 2001) and also that the radio morphology 
(typically classes FR I and FR II) may depend on the ab- 
solute magnitude of the host galaxy (Owen & Laing 1989; 
Owen & White 1991). This could be related to the differ- 
ent BH mass (proportional to the galaxy luminosity) and 
to the jet power (GhiscUini & Celotti 2001). 

As a whole, these studies have shown that the optical 
morphological/structural and photometrical properties of 
radio and non-radio ellipticals are remarkably similar, sug- 
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gesting that all ellipticals may go through a phase of nu- 
clear activity lasting for a small fraction of the total life 
of the galaxy. 

Few previous observations were devoted to investiga- 
tion of the dynamical properties of Radio Galaxies (RG). 
Apart a very small number (11) of radio sources included 
in the original work on the FP of ellipticals by Faber 
et. al. (1989) (hereafter FA89), only Smith, Heckman & 
Illingworth (1990; hereafter SHI90) performed a system- 
atic study on the stellar dynamics of Powerful Radio 
Galaxies (PRG). They compared </ie>, r^, and CTc mea- 
surements for a compilation of PRG with the distribution 
in the FP of the bright ellipticals studied by FA89, con- 
cluding that the FP of PRG is consistent with that of 
normal ellipticals. However they also found evidence for 
smaller than normal velocity dispersions and significant 
rotational support in galaxies with marked morphological 
peculiarity. 

In this paper we present a much deeper investigation 
of the FP of radio galaxies, based on a sample of 73 ra- 
dio galaxies. The data are in part new (22 objects) and in 
part collected from the literature (51 objects). Most of lit- 
erature values were derived from the Hypercat database 
(Prugniel & Maubon 2000). All values of </ie>, fe, and 
(7c have been processed in order to make them homoge- 
neous to a common standard (see Sec 4.2). Throughout 
the paper we assume i/o=50/i^^km s^^Mpc^^ and q^—Q. 



Table 1. Log of the Observations 
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2. Observations and data analysis 

We obtained medium resolution optical spectra of RG se- 
lected from the brightest {ran <15) objects in the sam- 
ple of 79 LzRG previously imaged, in the R filter, by us 
(Fasano et. al. 1996, Govoni et. al. 2000a, b). These ob- 
servations were aimed at deriving the velocity dispersion 
from stellar absorption lines. In Table ^ we give the list 
of the objects observed together with exposure times and 
position angles used for the observations. 

Optical spectroscopy was obtained in March/ April 
1998 and November 1998 with the ESO/Danish 1.52 m 
telescope at La Silla, using the Danish Faint Object 
Spectrograph and Camera (DFOSC). We used a CCD 
Loral/Lesser, with 2052 x 2052 pixels combined with an 
Echelle Grism of 316 grooves/mm yielding a velocity res- 
olution (FWHM) of 71 kms^i (for a sht 1" wide) in the 
range AA=4800-5800 A. A long slit, 2.0" wide and centered 
on the galaxy, was oriented along the apparent major axis 
of the radio galaxy. With this configuration we reach a ve- 
locity dispersion resolution Acr=60 kms~^; the scale per- 
pendicular to the dispersion is 0.39"/pixel. Template ref- 
erence spectra of standard stars of spectral type from G8- 
III to Kl-III, with low rotational velocity (Vxsin(i) <17 
kms"""^), were secured at the beginning and at the end 
of each night. During the observations the seeing ranged 
between 1" and 1.5". 

Optical spectra were reduced using standard proce- 
dures available in the IRAF package and includes bias 
subtraction, flat fielding, and wavelength calibration. The 



accuracy of the latter procedure was checked with mea- 
surements of the night sky Ao= 5577.32 A emission line. 
In all cases a precision of ^ 20 kms~^ was reached. In 
order to increase the signal to noise and to match the ob- 
served spatial resolution (assuming the mean seeing) with 
the plate scale across dispersion, the spectra were rebinned 
over three pixels perpendicular to the dispersion, obtain- 
ing an effective spatial resolution of 1.17"/pixel. The ob- 
served spectral range includes the Mg2 band (Ao= 5175.4 
A), the E-band (5269 A) and the Fel line (5335 A). 

The systemic velocity, corrected to the Sun, and 
the velocity dispersion a were determined using the 
Fourier Quotient method (Sargent et. al. 1977, Bertola 
et. al. 1984). The spectra were first normalized by sub- 
tracting the continuum, converted to a logarithmic scale 
and then multiplied to a cosine bell function that apodizes 
10% of the pixels at each end of the spectrum. This forces 
the ends of the spectra smoothly to zero. Finally the 
Fourier Transform of the galaxy spectra were divided by 
the Fourier Transform of a template star whose spectra 
(late G and early K spectral type) matched that of the 
galaxy. These spectra are used as templates of zero veloc- 
ity dispersion. 

The best-fit stellar template then yields a profile of the 
velocity dispersion a and the stellar velocity curve with 
relative errors for the galaxy. The r.m.s. of the determi- 
nations obtained with different template stars turned out 
to be less than 20 kms~^ for a and ~10 kms~^ for the 
systemic radial velocity Vr. The average values of these 
determinations were adopted as final values of a and Vr. 
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Fig. 1. Our velocity dispersion a measurements compared 
with the literature data 

Since early-type galaxies exhibit some gradients in the 
radial velocity and velocity dispersion (Davies et. al. 1983, 
Fisher et. al. 1995), the derived 'central' parameter CTc de- 
pends on the distances of the galaxies and the size of the 
aperture used for the observation. In order to compare 
our velocity dispersions with the data available in the lit- 
erature we applied aperture corrections according to the 
procedure given by J0rgensen et. al. (1995). The individ- 
ual measurements of a are therefore corrected to a circular 
aperture with a metric diameter of 1.19/i~^ kpc , equiva- 
lent to 3.4" at the distance of the Coma cluster. 

For 8 of the 22 RGs observed previous determina- 
tions of the velocity dispersion were found in the Hypercat 
database. The comparison between our and previous mea- 
surements of a is shown in Figure |^. When more than one 
measurement of a is present in the literature, we took the 
mean value quoted by Hypercat. The average difference 
between our and Hypercat values of a is: 

< o'our-o'Hyp > ~ 4±6.4 kms^^; r.m.s. — 18 kms^^ 

Additional notes on individual galaxies are given in the 
Appendix. 

Table |2| reports our measurements of CTc and the 
estimated uncertainities, together with the parameters 
</ie> and Te , taken from Govoni et. al. (2000a), that 
are relevant to construct the FP. The average surface 
brightness </ie> has been evaluated from the formula: 

</ie> = Rt + 51og(r^') -I- 2.51og(27r), 

where is the total R apparent magnitude and r" is the 
effective radius. All these quantities are corrected for the 
contribution of the point source (see Govoni et al. 2000b), 
for cosmological dimming and for i^T-correction accord- 
ing to Poggianti (1997). The distance/systemic velocities 
used to derive and </ie> are relative to the Cosmic 
Microwave Background (CMB) reference frame and come 



Table 2. Photometrical and kinematical Data for our 
sample of radio galaxies 
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from the Lyon Meudon Extragalactic Database (LEDA, 
Paturel et. al. 1997). 

3. Extended sample of LzRG 

We have collected from the literature stellar velocity dis- 
persion measurements of LzRGs for which photometric 
and structural parameters are also available. To make 
these data homogenous with our measurements we pro- 
cessed them applying the same procedure adopted for our 
subsample of LzRG (see above section and Govoni et. 
al. 2000b). 

We excluded from the analysis the radio galaxies be- 
longing to dumbbell systems or close galaxy pairs, as well 
as those having recession velocity (corrected to CMB) less 
than 3000 kms~^, since these circumstances may induce 
significant uncertainties in the parameters considered. 

3.1. Previous work on the Fundamental Plane of 
LzRGs 

An early investigation of the FP of LzRG was undertaken 
by FA89 in their study of early-type galaxies. They give 
</ie> and (in the B band), as well as a for 11 RGs. 
However 6 objects are dumbbell systems or have small 
recession velocity and for this reason are not considered 
in our analysis. For the remaining 5 we obtained <^e> 
in the R band applying a color correction based on the 
integrated color (B-R) obtained from Hypercat. 

The first systematic study of the FP of RGs was per- 
formed by Smith, Heckman & lUingworth (1990, SHI90) 
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who reported kinematical and photometrical data for 20 
powerful radio galaxies previously imaged in the V band 
(Smith & Heckman 1989, hereafter SH89), as well as for 

16 more galaxies from Heckman et. al. (1985). Among the 
latter ones, only 7 are supplied with complete photometric 
information, whereas 5 objects, out of the remaining 27 in 
the total sample, turn out to belong to dumbbell systems 
or close pairs. This reduces the usable sample of SH190 
to 22 objects available for FP analysis. In order to derive 
</ie> consistently with our definitions (see above formula) 
we used their isophotal magnitudes V25 and the effective 
radii reported in SH89 (Table 8) and applied a correc- 
tion of -0.16 (average difference between mjot and m24.5 
in our study of RG; Govoni et. al. 2000b) to obtain total 
magnitudes. Then we corrected for the color term (V-R) 
using the colors from SH89, when available, or assuming 
(V-R)=0.55 in the other cases. 

3.2. Additional data for the Fundamental Plane of 
LzRG 

In order to improve the statistics of our analyses, we have 
searched the literature for additional LzRG that have mea- 
surements of velocity dispersion, <^e> and rg. We have 
combined photometric data of LzRGs from the studies 
of Ledlow & Owen (1995, L095) and Gonzalez-Serrano 
& Carballo (2000, GCOO) whith the velocity dispersions 
given by Wegner et. al. (1999, EFAR) or by Hypercat. 
The cross identification betwen L095 and EFAR added 

17 new LzRGs to our compilation, whereas that between 
GCOO and Hypercat adds 7 more LzRGs. In both cases, in 
order to derive </ie> consistently with our definitions, we 
used the magnitudes 77124.5 and the effective radii re given 
in the original papers and applied the same corrections as 
for the data of SHI90. 

In total, the sample of LzRG for which we collected 
data from the literature consists of 51 objects. The rel- 
evant data for these objects are given in Tables || and 
^. These, together with the sample for which we present 
new data, leads to a total of 73 LzRG, and represents the 
largest dataset of photometric and spectroscopic data on 
Radio Galaxies till now. 



3.3. Comparison sample of non-radio elliptical galaxies 

In order to compare the FP of LzRGs with that of nor- 
mal galaxies, we used the sample of radio-quiet early-type 
galaxies studied by J0rgensen et. al. (1996, JFK96). It 
is by far the richest, still homogeneous sample for which 
both photometric and kinematic information is available. 
Previous samples of early type galaxies used to describe 
the FP (e.g. FA89) have a significantly larger scatter with 
respect to the data of JFK96. Since the JFK96 sample 
consists of cluster galaxies, it has also the advantage of 
the reliability of the distance determination. 

To compare our data on LzRG with the JFK96 data, 
we converted their Gunn r photometry to the R band, 



Table 3. Radio Galaxies from literature 
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by means of the average relation R — r=0.3, given in 
J0rgensen (1994). 

4. The Fundamental Plane of radiogalaxies 

In order to derive the parameters describing the FP: 

log r'^P" = a log o- + (3<fie> r-1, (1) 

we minimized the root square of the residuals perpendic- 
ular to the plane. This procedure is to be preferred with 
respect to the classical one (minimizing the root square of 
the residuals along a given axis) when there are measure- 
ment errors in all observed quantities. It is also slightly 
different from that used by JFK96, which minimizes the 
sum of the absolute residuals perpendicular to the plane. 

In Table || the FP coefficients a, (3 and 7 obtained us- 
ing this fitting procedure for different subsamples of radio 
and non radio galaxies are reported. Since the adopted 
fitting procedure does not provide analytical form for the 
uncertainties of the parameters, we indicate in Table 
the la uncertainties computed according to the classical 
formalism and assuming the values of the FP coefficients 
given in the table. 

Each sample in the table is indicated with a capital 
letter and the fits obtained merging together two or more 
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Table 4. Radio Galaxies from literature, continue 
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samples are indicated "with sequences of letters, according 
to the previous correspondence (for instance JO means 
JFK96+This Work). 

We note from Table || that, in spite of the slightly dif- 
ferent fitting procedure, the FP coefficients we found for 
the JFK96 sample of radio-quiet galaxies (our compari- 
son sample) are indistinguishable from those given in the 
original JFK96 paper (see their equation [1]). In Figure || 
we plot all the galaxies in the above described samples, 
together -with our FP fit of non radio ellipticals (JFK96). 

From Tabic || the FP coefficients derived from our sam- 
ple of 22 radio-galaxies turn out to be consistent (-within 
Icr of the estimated uncertainties) -with those relative to 
the control sample of non radio galaxies. The same hap- 
pens for the 17 radio-galaxies in the L095 sample. On 
the contrary, fitting the SHI90 sample of 22 radio-galaxies 
yields FP coefficients different (-well beyond Scr) from those 
relative to the comparison sample. This is likely to be 
due to the presence of a fe-w outliers in their galaxy sam- 
ple, as sho-wn in Figure |^. Finally, given the small size of 
FA89 and GCOO samples a reliable estimates of the FP 
can not be derived from them. Nevertheless, we have re- 
ported in Table || the coefficients resulting from the for- 
mal fit of these samples. In the case of the GCOO sample, 
they are consistent -with those relative to the control sam- 
ple, -whereas considerable differences (again caused by the 
presence of an outlier) come in the case of the FA89 sam- 
ple. 



In Table g| we also report the fits obtained merging to- 
gether f-wo or more galaxy samples. Among them, -we as- 
sume as representative of the global FP of elliptical galax- 
ies, the fit: 

log^fcpc ^ 1.27 log cr + 0.326<Aie>fl - 8.56 (2) 

-which includes the comparison sample (JFK96) and the 
radio samples from this -work(O), from L095(L) and from 
GCOO(G). This fit allows to extend the FP up to the bright 
end of the luminosity function of early-type galaxies. 

The M/L ratio is a function of the stellar population 
and dark matter content of galaxies. In Figure ^ we show 
the relation between the Mass to Light ratio (M/L) and 
the velocity dispersion a. We derived the masses of our 
galaxies by using the relation M=5reCr^/G (Bender et. 
al. 1992). This figure indicates that LzRGs follow the same 
relation found for normal elliptical galaxies. The scatter 
of the M/L ratio as a function of log a is ■~0.2dex, similar 
to that found by JFK96. 

5. Conclusions 

The main conclusion of this study is that the fundamen- 
tal plane of radio galaxies, as defined by our collection of 
data for 73 objects, is consistent with the one defined by 
normal, non-radio elliptical galaxies. Some recent results 
on BH demography (Ferrarese & Merrit 2000, Gebhardt 
et. al. 2000) strongly suggest that all galaxies may host a 
massive BH in the nucleus, and that the BH mass is pro- 
portional to the mass of the spheroidal component. This 
means that, virtually, all ellipticals have the basic ingre- 
dient for becoming active. Considering the small amount 
of gas (few solar masses per year) necessary to sustain the 
radio emission of even the most powerful radio sources, the 
activity could be triggered by extremely modest alteration 
of the status of equilibrium in which galaxies settled soon 
after they form. According to this view it is therefore not 
surprising that radio and non radio ellipticals have indis- 
tinguishable global properties, irrespective of their nuclear 
activity. Here, we present one more result supporting this 
scenario. 

Both metallicity effects and age variations have been 
suggested as the main cause (Faber et. al. 1995) of the 
strong dependence of M/L on cr, that is on the total mass 
of the galaxy. This relation is consistent with a progres- 
sive reddening of the stellar population going from small 
to big galaxies (e.g., Prugniel & Simien 1996), but most of 
dependence remains unexplained. It looks like that ellipti- 
cals pass from being baryon dominated to be dark matter 
dominated with increasing luminosity, with a "fine tun- 
ing" required by the small and basically constant disper- 
sion around the FP. 

Appendix A: Notes on individual galaxies 

0055-016 - UGC 595 - 3C29 this galaxy belongs to the rich 
cluster A119 and is D26 in the Dressier (1980) list. Smith 
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Fig. 2. The Fundamental Plane of various samples of LzRG (see legend and Table |5|) and normal (non radio) ellipticals 
(JFK96), compared with the best fit (dotted line) to the JFK96 data 



Table 5. The FP in R 



Sample 


Ngal 


a 




7 


rms 


J (JFK96) 


229 


1.24±0.04 


0.330±0.007 


8.58±0.06 


0.059 


(This Work) 


22 


1.36±0.23 


0.295±0.031 


8.15±0.23 


0.071 


L (L095) 


17 


1.24±0.23 


0.311±0.015 


8.21±0.18 


0.047 


G (GCOO) 


7 


2.55±0.74 


0.311±0.073 


11.37±0.53 


0.055 


S (SHI90) 


22 


2.60±0.31 


0.161±0.021 


8.19±0.26 


0.070 


F (FA89) 


5 


2.31±0.48 


0.072±0.043 


5.86±0.50 


0.030 


OLG 


46 


1.58±0.16 


0.310± 0.015 


8.99±0.15 


0.063 


OLGSF 


74 


1.92±0.15 


0.256±0.013 


8.69±0.15 


0.082 


JOLG 


275 


1.27±0.04 


0.326±0.007 


8.56±0.06 


0.060 


JOLGSF 


303 


1.35±0.04 


0.308±0.007 


8.40±0.06 


0.069 



et. al. (1990) measured a velocity dispersion of 199±18 
kms~^but recently Wegner et. al. (1999) give a value of 
253±17 kms^^ in better agreement with our measure- 
ments. 



0131-360 - NGC 612 This galaxy show a strong dust 
lane along the apparent major axis, our spectrum has been 
obtained perpendicularly to the dust lane, i.e. along the 
minor axis of the optical galaxy, no rotation is visible. 
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Fig. 3. The M/L ratio versus the velocity dispersion relationship of Radio Galaxies compared with that for early type 
galaxies (JFK96) 



0915-118 - 3C218, Hydra A, A well studied galaxy 
showing strong [OIII](A = 5007 A) emission lines. 
Heckman et. al. (1985) found a = 308 ±38 kms~^, slightly 
higher than our measurement 278 ± 78 kms~^. 

1107-372 - NGC 3557, We measure a velocity disper- 
sion of 295±7 kms~^ and a velocity gradient of '--^150 
kms~^ in the central 5" in good agreement with data 
available in Hypercat. 

1258-321 - ESO443-G024, this galaxy is in the cluster 
A3537. Dressier et. aZ.(1991) measured a velocity disper- 
sion of 279±27 kms~^ in good agreement with our data. 

1333-337 - IC 4296, Our value of sigma is lower 
than the mean value quoted by Hypercat (340 kms^^). 
The agreement is better with the measure of Franx et. 
al. (1989). From our spectrum, taken at PA=60° we mea- 
sure a rotation of ^-^71 kms^^ in agreement with data 
taken at the same position angle by Saglia et. al. (1993). 



1400-337 - NGC 5419, Our measure is in good agree- 
ment with the mean value quoted by Hypercat. The spec- 
trum, taken at PA=84°, show a maximum rotational ve- 
locity of ~90 kms~^ over 5" . 

1514-h072 - UGC9799, our measure of the velocity dis- 
persion is very close to the mean value of 247 kms~^ 
quoted in Hypercat. 
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